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ABSTRACT 


The  problem  of  intersecting  N half-spaces  in  K space  is  transformed  to-SaK-*- 
. problem^  of  constructing  the  convex  hull  of  N points  in  K space  and  a simple 
intersection  problem.  This  enables  one  to  intersect  the  N K-dimensional  half  spates  in 
0(K*H(N,K)J  time,  wht^e  H(N,K)  is  fhe  time  required  to  construct  the  convex  hull  of  N 
points  in  K space.  Fo^  two  and  three  dimensions  the  algorithm  takes  OfNIogN)  time  in 
the  worst  case,  but  untfpr  fairly  robust  conditions  the  expected  time  is  only  O(N).  It  is 
also  shown  that  an  algorithm  for  intersection  of  half  spaces  can  be  used  to  construct 
the  convex  hull  of  points  in  K space.  Thus,  the  intersection  of  half  spaces  and  convex 
hull  of  points  problems  are  essentially  equivalent. 


This  research  was  partially  supported  by  the  Office  of  Naval  Research  under  contract 
number  N00014-76-C-0829. 
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1.  Introduction 

There  has  been  much  recent  work  on  the  problem  of  intersecting  half 
spaces.  Shamos[75]  and  Shamos  and  Hoey[76]  describe  an  CXNIogN)  time  algorithm  in 
two  dimensions,  and  Zolnowsky{77]  and  Preparata  and  Muller[77]  have  developed 
OfNIogN)  time  algorithms  for  three  dimensions.  These  algorithms  are  optimal  within  a 
constant  factor  because  there  is  an  fl(NlogN)  lower  bound  for  two  or  more  dimensions 
(Shamos[75]>. 

We  will  show  that  an  algorithm  for  constructing  the  convex  hull  of  a set  of 
points  in  K space  can  be  used  lo  intersect  half  spaces  in  K space,  and  vice-versa.  In 
fact,  within  a factor  of  K,  the  upper  and  lower  bounds  for  one  problem  apply  to  the 
other.  Thus,  (much  of)  the  knowledge  concerning  convex  hulls  of  points  in  K space 
applies  to  the  intersection  of  half  spaces  in  K space.  Graham[72]  and  Preparata  and 
Hong[77]  describe  CXNlogN)  time  algorithms  for  constructing  the  convex  hull  of  N points 
in  the  plane.  Bentley  and  Shamos[78]  and  Eddy[77]  give  fast  expected  time  algorithms 
in  the  plane.  Preparata  and  Hong[77]  have  also  produced  an  0(NlogN)  time  algorithm 
for  the  convex  hull  of  N points  in  3 space.  However,  in  four  dimensions  there  is  an 
O(N^)  lower  bound  because  the  convex  hull  can  have  0(N^)  edges  (Grunbaum[67J, 
p.193).  Chand  and  Kapur(70]  describe  a convex  hull  algorithm  for  an  arbitrary  number 
of  dimensions  whose  complexity  has  yet  to  be  analyzed. 

This  paper  will  first  describe  the  two  dimensional  problem  and  the  solution 
obtained  by  use  of  the  transform.  Then  we  move  up  to  three  dimensions,  using  the 
transform  to  solve  the  problem  in  O(NlogN)  time.  Finally,  four  (and  more)  dimensions 
'will  be  tackled,  followed  by  some  fast  expected  time  algorithms  and  applications.  All 
. algorithms  in  this  paper  use  the  RAM  (with  reals)  model  of  computation. 


2.  The  two  dimensional  problem 

Shamos[75]  and  Shamos  and  Hoey[76]  show  that  the  intersection  of  N half 
planes  has  time  complexity  0(NlogN).  Their  algorithm  for  constructing  the  intersection 
**  in  OfNIogN)  time  relies  on  Shamos'  CXN)  time  algorithm  for  intersection  of  convex 
polygons.  In  this  paper  an  entirely  different  0(NlogN)  time  algorithm  is  described. 

I 

i 

In  Figure  1 the  intersection  of  N half  planes  is  indicated  by  the  shaded 
region.  The  half  planes  are  partitioned  into  two  sets,  UPPER  and  LOWER.  A half  plane 
is  in  set  UPPER  if  the  line  at  its  boundary  is  above  the  rest  of  the  half  plane.  Similarly, 
a half  plane  is  in  the  set  LOWER  if  Ihe  line  at  its  boundary  is  below  the  rest  of  the  half 
plane.  (If  any  boundary  lines  arc  vertical,  then  rotate  all  N half  planes  a small  angle.) 
The  reason  for  producing  this  partition  will  not  be  explained  now,  but  it  will  be 
apparent  later.  Thus,  the  problem  is  now  divided  into  three  parts:  (1)  the  intersection 
of  the  UPPER  half  planes,  (2)  the  intersection  of  the  LOWER  half  planes,  and  (3)  the 
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Figure  1:  Intersection  ot  N half  planes.,  Figure  2:  Intersection  of  regions  U and  L. 


As  shown  in  Figure  2,  part  (3)  is  relatively  easy.  If  U and  L in  Figure  2 have 
0(N)  vertices  then  the  intersection  (shaded  region)  can  be  constructed  in  0(N)  time. 
The  algorithm  is  described  in  detail  in  ALGORITHM  1NTERSECTCHAINS  below. 


ALGORITHM  INTERSECTCMAINS 

Input:  Integers  Nl,  N2  (where  N - N1  ♦ N2)  and  real  vectors  UX[0:N1*1],  Uy[0:Nl*l], 
LX[0:N2+1  ],  L [0:N2*1]  such  that 

-<»  - UxfO?  < Uxfl]  <...<  UxfNl  1 < UJNl  *1  ]-  ® 

-oo  - LX10]  < Lxll]  < . . . < LXIN2]  < LX|N2  + I]  - oo. 

Output:  Integer  K (number  of  vertices  on  the  intersection),  HX[1:K],  FL[1:K). 

Time:  0(N),  Space:  0(N).  y 

(1)  Scan  U and  L (vectors  Ux,  Uy,  Lx,  and  Ly)  from  left  to  right  until  two  segments 
intersect  at  a point  P.  (If  no  segments  intersect  then  the  intersection  of  U and  L is 
empty.)  The  scan  can  bo  done  in  CXN)  time  in  a manner  similar  to  the  0(N)  time  merge 
in  the  merge  sort  algorithm. 

(2)  Scan  U and  L from  right  to  left  until  two  segments  intersect  at  a point  Q. 

(3)  If  P »*  0,  then  return  (in  vectors  and  Hy)  the  concatenation  of  the  chains  of  line 
segments  of  U and  L between  points  P and  Q. 

(4)  If  P - Q,  then  the  intersection  is  unbounded  (or  just  the  point  P - Q).  In  the  case 
of  an  unbounded  intersection  one  must  determine  whether  to  return  the  chains  to  the 
left  of  P or  the  chains  to  fhe  right  of  P.  This  can  be  determined  by  comparing  tire 
slopes  of  the  rays  bounding  the  left  and  right  sides  of  U and  L.  If  the  slope  of  the  left 
ray  of  U is  less  than  the  slope  of  the  left  ray  of  L,  then  return  the  chains  to  the  left  of 
P.  Otherwise,  return  the  chains  to  the  right  of  P. 
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This  leaves  parts  (I)  and  (2),  the  intersection  of  the  UPPER  half  planes  and 
the  intersection  of  the  LOWER  half  planes.  Only  (1)  will  be  described  since  (2)  can  be 
solved  similarly.  Assume  that  the  N half  planes  are  all  UPPER  half  planes.  Some  of 
these  half  planes,  such  as  half  plane  k in  Figures  3a  and  3b,  do  not  bound  any  side  of 
the  region  o<  intersection.  It  would  be  nice  to  be  able  to  find  all  such  half  planes  and 
throw  them  away  since  they  do  not  contribute  to  the  final  result.  Once  that  is  done 
one  can  find  the  intersection  of  the  UPPER  half  planes  rather  easily.  As  one  can  see  in 
Figure  2,  the  slopes  of  the  sides  of  the  chain  U are  rnonolonic  decreasing  as  one 
travels  from  left  to  right.  Thus,  given  the  lines  determined  by  the  sides,  one  need  only 
sort  the  lines  by  slope  lo  determine  the  order  in  which  they  intersect  to  form  Ihe 
sides.  The  sort  costs  CKNIogN)  time  so  once  all  the  redundant  half  planes,  such  as  k in 
Figures  3a  and  3b,  are  removed,  Ihe  intersection  of  Ihe  UPPER  half  planes  can  be 
determined  in  0(NlogN)  lime. 


Figure  3:  (a)  & (b)  - k is  redundant,  (c)  - k is  nonredundant. 


How  does  one  determine  which  half  planes  are  redundant  and  which  are  not? 
There  are  two  conditions  which  need  to  be  checked. 

(1)  Line  k (Figure  3)  must  be  above  the  point  P where  lines  i and  j meet. 

(2)  The  slope  of  line  k must  lie  between  the  slopes  of  lines  i and  j.  That  is, 
slope(i)  < slope(k)  < slope(j)  or  slope(j)  < slope(k)  < slope(i).  This  is  true  if  i and 
j have  slopes  of  different  sign  (Figure  3a)  or  ol  the  same  sign  (Figure  3b). 

Figure  3c  shows  what  will  happen  if  the  first  condition  is  satisfied  but  the  second  is 
not.  A half  plane  k is  therefore  redundant  iff  there  exist  half  planes  i and  j such  that 
the  two  above  conditions  arc  satisfied. 

How  fast  can  one  determine  (non)redundancy  for  each  of  Ihe  N UPPER  half 
planes?  Certainly  one  approach  is  to  test  all  pairs  of  half  planes  i and  j for  each 
helfplane  k.  That  costs  CHN^)  lime,  though,  which  is  not  good.  One  could,  perhaps, 
modify  this  approach  and  obtain  a (aster  algorithm.  However,  there  is  an  enlirely 
different  way  of  looking  at  this  problem  that  is  clearly  better  and  that  is  the  approach 
which  will  be  taken. 
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3.  The  Transform 

The  transform  exploits  a natural  duality  between  points  and  lines  in  the 
plane.  For  an  arbitrary  point  (x,y),  consider  the  set  of  all  lines  in  slope-intercept  form 
which  pass  through  the  point.  (Thus  this  won’t  work  for  vertical  lines.)  If  a line  is 
represented  by  an  ordered  pair  of  the  form  (slope, intercept),  then  this  set  is 
{ (a,b)  | y ■ a**  ♦ b ).  Note  that  this  transforms  a point  (x,y)  to  a line  b ■ -x*a  * y. 
Note  also  that  the  representation  of  a line  y - a*x  b in  the  form  (a,b)  is  a natural 
transform  of  a line  to  a point.  Thus,  points  transform  to  lines  and  lines  transform  to 
points  by  the  formulas 

y - a*x  ♦ b -*  (a,b),  and 
(x.y)  -*  b - -xta  ♦ y. 

This  transform  has  an  interesting  property:  Distances  in  the  y-c oordmate 
between  points  and  lines  are  preserved.^  The  difference  in  y coordinate  between 
point  (c,d)  and  line  y • e*x  ♦ f is  d - (e*c  ♦ f).  The  difference  in  the  transforms  b - - 
c*a  + d and  (e,f)  is  (-c*e  ♦ d)  - f,  which  is  the  same.  It  follows  from  this  that  incidence 
is  preserved.  If  point  (c,d)  is  on  line  y - e*x  ♦ f,  then  it  holds  also  for  their 
transforms  - point  (e,f)  is  on  line  b ■ -c*a  ♦ d.  Note  further  that  not  only  is  the 
magnitude  of  the  distance  (in  the  y-coordinate)  preserved  but  also  its  sign.  Thus, 
above/below -ness  is  preserved.  If  (c,d)  is  above  (below)  line  y - e*x  + f,  then  the 
transform  of  (c,d)  is  above  (below)  the  transform  of  y - e*x  ♦ f. 

There  is  another  property  of  the  transform  which  should  be  mentioned.  The 
transform  is  not  involutory,  but  composition  of  it  four  times  produces  the  following: 

(x,y)  -*  b - -x*a  ♦ y -»  (-x,y)  -♦  b - x*a  ♦ y -»  (x,y) 

Only  a slight  change  is  required  to  make  the  transform  its  own  inverse.  Express  lines 
in  the  form  y ♦ a*x  ♦ b ■ 0 rather  than  y • a*x  ♦ b.  Then  it  is  true  that  y + a*x  ♦ b ■ 
0 ♦♦  (a,b).  But  this  has  the  unfortunate  side  effect  that  above/below-ness  between 
points  and  lines  is  not  preserved;  it  is  reversed.  If  point  (c,d)  is  above  line  y ♦ e*x  ♦ f 
■ 0 then  the  transform  of  (c,d)  will  be  below  the  transform  of  line  y ♦ e*x  ♦ f - 0. 


*The  restriction  in  the  y coordinate  is  important  because  it  can  be  shown  to 
be  impossible  to  preserve  the  Euclidean  distance  between  a point  and  a tine  under  a 
duality  transform. 

There  are  other  duality  transforms  which  preserve  incidence,  such  as 
Plucker’s  transform  [Shamos  77J. 
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4.  Application  of  the  transform  to  the  two  dimensional  problem 

Now  it  will  be  shown  how  the  transform  enables  one  to  intersect  the  UPPCR 
(or  LOWER)  halt  planes  Iasi.  More  specifically,  the  transform  enables  an  efficient 
mechanism  for  eliminating  the  “redundant*  half  planes.  Recall  the  two  conditions  for 
redundancy  of  an  UPPER  half  plane:  » half  plane  K is  redundant  iff  there  exist  half 
planes  i and  j such  that  (H  line  k is  above  the  point  P where  lines  I and  j intersect, 
and  (2)  the  slope  of  k is  between  the  slopes  of  lines  i and  j.  In  the  ab  plane  there  is  a 
corresponding  interprelahon. 


r *k.i 


Figure  A:  Transform  of  a redundant  half  plane. 


In  Figure  4,  line  k is  above  point  P in  the  xy  plane.  This  is  transformed  to  a 
point  k which  is  above  line  P in  the  ab  plane.  (Above /below -ness  between  lines  and 
points  is  preserved  by  the  transform.)  The  slope  of  a line  in  the  xy  plane  is  the  a 
coordinate  of  the  corresponding  point  in  the  ab  plane.  Thus,  as  in  Figure  A,  line  k with 
a slope  between  the  slopes  of  lines  i and  j transforms  to  a point  k with  a coordinate 
between  the  a coordinates  of  points  i and  j in  the  ab  plane.  Figure  5 shows  the  result 
of  applying  this  procedure  to  several  half  planes.  In  the  ab  plane  the  points 
connected  by  line  segments  are  the  transforms  of  the  "nonredundant"  half  planes.  The 
points  above  these  segments  correspond  to  the  redundant  half  planes. 


Figure  5:  Transform  of  N UPPER  half  planes. 


Look  closely  at  those  line  segments  in  the  ab  plane.  They  are  exactly  the 
bottom  part  of  the  convex  hall  of  the  points  in  the  ab  plane!  The  problem  of 
intersecting  N UPPER  half  planes  has  been  reduced  to  the  problem  of  constructing  the 
(bottom  part  of  the)  convex  hull  of  N points.  The  convex  hull  of  N points  in  the  plane 
can  be  constructed  in  CXNIogN)  time  (Graham[72]).  This  leaves  only  the  detail  of 
separating  the  top  from  the  bottom  part  of  the  hull.  To  do  that,  find  the  leftmost  point 
of  the  hull  in  0(N)  time.  Then  traverse  the  hull  on  the  bottom  side  until  the  rightmost 
point  is  reached. 
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The  problem  of  internetting  N half  planes  has  been  broken  into  three  parts. 
Part  (1),  the  intersection  ol  the  UPPER  half  planes,  has  just  been  shown  to  cost  only 
O(NlogN)  time.  Part  (2),  the  intersection  of  the  LOWER  half  planes,  is  equivalent  to  part 
(1)  so  it  can  also  be  done  in  (XNIogN)  time.  Part  (3),  the  intersection  of  the  results  of 
parts  (1)  and  (2),  has  been  shown  to  cost  only  0(N)  lime.  Thus,  the  entire  problem  tan 
be  solved  in  (XNIogN)  time. 

It  is  interesting  to  note  that  one  can  also  use  an  intersection  of  half  planes 
algorithm  to  produce  a convex  hull  ol  points  algorithm.  First  transform  all  N points  to 
UPPER  half  planes  by  the  lormula 

(x,y)  -»  b - xta  ♦ y 

and  intersect  the  half  planes.  Transform  back  to  obtain  the  lower  part  of  the  hull. 
Then  transform  all  N points  to  LOWER  half  planes  and  intersect  the  half  planes. 
Transform  back  for  the  upper  part  of  the  hull.  Merging  the  upper  and  lower  parts  is 
trivial,  since  the  leftmost  and  rightmost  points  will  be  in  each  one.  Total  time  is 
CXNtogN). 


5.  The  three  dimensional  problem 

The  technique  just  used  in  two  dimensions  can  be  extended  to  three 
dimensions.  The  half  spaces  are,  as  before,  first  partitioned  into  the  two  sets  UPPER 
and  LOWER.  A half  space  is  in  the  UPPER  :et  if  the  plane  at  its  boundary  is  above  the 
rest  of  the  half  space.  Similarly  lor  LOWER.  The  problem  ol  intersecting  the  half 
spaces  is  then  broken  into  the  subproblems  (1)  intersect  the  UPPER  half  spaces,  (2) 
intersect  the  LOWER  half  spaces,  and  (3)  intersect  the  results  of  (1)  and  (2).  This 
schema  will  still  work  in  three  dimensions,  but  the  algorithm  for  part  (3)  is  more 
complicated  and  does  not  seem  to  generalize  well  to  higher  dimensions.  To  avoid  these 
complications  a slightly  different  method  for  intersecting  half  spaces  in  three 
dimensions  will  be  presented. 

Recall  that  the  partition  of  the  half  spaces  into  the  sets  UPPER  and  LOWER  is 
determined  by  a comparison  in  the  z coordinate.  The  choice  of  the  z coordinate  is 
totally  arbitrary.  Any  of  the  3 coordinates  will  do  just  as  well. ^ Thus,  there  are  3 
possible  sets  UPPER  and  LOWER.  To  distinguish  them  from  each  other,  use  the  notation 
UPPER j and  LOWER j for  the  first  (x)  coordinate,  UPPER2  and  LOWER2  for  the  second 
(y)  coordinate,  and  UPPER3  and  LOWER3  for  the  third  (z)  coordinate.  Also,  U and  L,  the 
intersections  of  the  UPPER  and  LOWER  sets,  will  now  be  written  Uj  and  lj  for  the  it h 
coordinate.  The  new  method  for  solving  the  intersection  of  U and  L (that  is,  U3  and 
L3),  uses  all  3 of  the  pairs  U;  and  Lj.  (See  Figure  6 for  an  illustration  of  the  three  sets 
U and  L for  an  octahedron.)  It  will  be  shown  how  the  extra  Uj  and  Lj  make  the 
Intersection  much  simpler.  Then  it  will  be  shown  how  to  construct  U3,  L3,  U2,  L2.  Uj, 
and  Lj. 


*ln  fact,  so  will  any  linear  combination  of  the  coordinates. 
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Figure  6:  The  three  sets  U and  L for  an  octahedron. 


An  example  of  the  two  dimensional  case  is  shown  in  Figure  7.  The  curves 
U2  and  1-2  of  Figure  7a  meet  at  two  points,  P and  0-  Once  P and  Q are  found,  it  is 
easy  to  construct  the  entire  polygon  defining  the  intersection  of  regions  U2  and  I2. 
(See  the  last  part  of  INTFWSECTCHAINS  in  Section  2.)  The  leftmost  intersection  point  P 
of  U2  and  L2  in  7(a)  is  simply  the  leftmost  point  of  Lj  in  7(b).  Thus,  rather  than  using 
algorithm  1NTERSECTCHAINS  (Section  2)  to  find  point  P,  one  can  simply  construct  Lj 
and  find  the  leftmost  point.  (Similarly  for  point  Q.)  Does  this  always  work?  No,  but  the 
exceptions  are  not  a serious  problem. 


Figure  7:  An  easy  way  to  find  P and  Q. 


Suppose  that  point  P is  not  only  the  leftmost  point  of  the  shaded  region,  hut 
also  the  highest  (or  lowest)  point  (Figure  8).  Then  P will  be  one  of  the  two 
intersection  points  for  not  only  U2  and  L2,  but  also  for  Uj  and  It.  Such  points  will  be 
called  extreme  intersection  points.  (The  cases  of  unbounded  or  null  intersections  will 
be  treated  later.)  In  fact,  as  shown  in  Figure  8,  P and  Q can  be  extreme  intersection 
points  simultaneously.  Thus,  in  the  two  dimensional  worst  case,  the  use  of  U}  and  Lj 
to  find  P and  Q does  not  seem  entirely  successful,  since  the  two  extreme  intersection 
points  may  still  be  left  unknown.  However,  in  higher  dimensions  the  situation  looks 
relatively  better. 


(a)  (b) 


Figure  8:  P and  Q are  extreme  intersection  points. 

Thm:  In  K dimensions  there  are  at  most  two  extreme  intersection  points  o(  N halt 
spaces. 

PI:  To  be  an  extreme  intersection  point,  the  vertex  must  be  extreme  in  all  K 

coordinates.  (That  is,  greatest  or  smallest  of  all  the  vertices  in  each  coordinate.)  It  is 

assumed  that  no  vertices  have  any  coordinates  in  common,  since  the  probability  of  that 

happening  is  zero.  Assume  that  P,  Q,  and  R are  distinct  extreme  intersection  points. 

Since  P and  Q are  extreme  intersections,  then  if  the  ith  coordinate  of  P is  the  smallest 

(largest)  over  the  vertices,  then  the  ith  coordinate  of  Q is  the  largest  (smallest). 

However,  the  ith  coordinate  of  R must  also  be  the  largest  or  smallest  i coordinate  and 

that  can’t  be  because  points  P and  Q exhaust  both  of  those  possibilities.  Thus  there 

i are  at  most  two  extreme  vertices.  QED. 

« 
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The  fact  that  there  are  at  most  two  extreme  intersection  points  is  used  in 
the  algorithm  (below)  for  intersecting  half  spaces  in  three  space.  This  algorithm 
assumes  that  one  can  construct  Uj,  lj,  Uj,  I2.  U3,  and  1 3 in  CXNlogN)  time.  That  will 
be  proven  later  in  this  paper,  but  lor  now  simply  assume  that  such  an  algorithm  exists 
and  It  will  be  shown  how  to  construct  the  intersection  of  U3  and  I3  (and  thus  the 
intersection  of  all  N half  spaces)  in  0<N)  additional  time. 

ALGORITHM  TO  INTERSECT  HALF  SPACES  IN  THREE  DIMENSIONS 

Input:  Integer  N > 0,  real  vectors  X[l:N],  Y[1:N],  Z[1:N], 

Output:  Integers  K (number  of  vertices  of  the  polytope),  E (number  of  edges),  F 
(number  of  faces),  real  vectors  H^liK],  R,[1:K],  HZ[1:K],  integer  arrays 
EndPointsOfEdges[  I :E,1 :21  FacesOfEdges[l  £,1:2],  EdgesOfFacesPointer[  I :F], 
EdgesOfFaces[  1 :2*E]. 

Time:  0(N!ogN),  Space:  0<N). 

(1)  Construct  U3,  I3,  U2,  13.  Uj.  and  Lj  in  CXNlogN)  time.  (Section  6 of  this  paper.) 

(These  structures  are  represented  the  same  as  the  output  for  this  algorithm.) 

(2)  For  the  most  part,  the  edges  where  U3  and  L3  meet  have  already  been  constructed 

as  edges  of  U2  and  13  (and  Uj  and  Lj).  For  each  of  the  0(N)  faces  of  and  L2 
one  can  determine  in  CK1)  time  whether  it  belongs  in  UPPER3  or  LOWER3.  Thus, 


mmmm 
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(a)  (b) 

Figure  9:  Dashes  indicate  boundaries  of  (a)  U3  and  L3  and  (b)  U2  and  1 2 

in  0(N)  time  one  ran  search  all  of  the  faces  of  U2  and  L2  and  determine  where 
the  UPPER3  faces  of  U2  and  l meel  the  IOWER3  faces.  The  edges  where  these 
faces  meet  are  edges  of  the  intersection  of  U3  and  I3.  (See  Figure  9.) 

(3)  Compare  U3  and  I3  with  Uj  and  Lj  as  in  (2)  in  CXN)  time.  By  the  above  theorem, 

there  will  be  at  most  two  vertices  of  U3  n L3  which  have  not  yet  been 
determined  in  steps  (2)  and  (3). 

(4)  If  the  chain  of  edges  connecting  U3  and  I3  (constructed  in  steps  (2)  and  (3))  is  a 

closed  loop  then  there  are  no  extreme  vertices  and  the  intersection  of  U3  and  I3 
is  complete. 

(5)  If  the  chain  of  edges  connecting  U3  and  I3  is  open  at  one  end  then  either  (a)  the 

intersection  of  U3  and  I3  is  an  unbounded  region,  or  (b)  there  is  one  extreme 
intersection.  If  the  chain  is  broken  into  two  parts  then  either  (c)  there  are  two 
extreme  intersections,  or  (d)  there  is  one  extreme  intersection  and  the 
intersection  of  U3  and  I3  is  unbounded.  For  these  cases  of  one  or  two  breaks  in 
the  chain,  apply  step  (6)  to  find  any  remaining  vertices. 

(6)  To  distinguish  an  extreme  intersection  point  from  an  unbounded  intersection:  (a) 

find  the  two  edges  at  the  two  ends  of  the  chain  and  call  them  A and  B,  (b)  find 
the  faces  above  and  the  faces  below  A and  B which  have  A or  B as  edges.  Call 
them  Aut  Aj,  By,  and  Bj.  (c)  Find  the  (infinite)  faces  at  the  intersection  of  Au  and 
Aj,  and  Bu  and  Bj.  Call  them  Aj  and  Bj.  (d)  In  CXN)  time  find  the  intersection,  if 
any,  of  Aj  and  Bj 

It  has  just  been  shown  that  given  U3,  L3,  U2,  L£,  Uj,  and  lj,  the  intersection 
of  Ug  and  L3  can  be  constructed  in  CXN)  time.  Now  it  must  be  shown  how  to  construct 
U3,  L3,  Ug.  t-2»  ^l*  Jn<1  in  0<MogN)  t'me.  Since  constructing  one  of  these  six 
intersections  is  essentially  the  same  as  constructing  any  other,  the  following  text  will 
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describe  just  the  intersection  of  a set  of  half  spaces  called  UPPER  and  the  intersection 
will  be  called  U. 

The  algorithm  for  constructing  U in  three  dimensions  is  analogous  to  the  t wo 
dimensional  case.  In  brief,  this  is  how  it  will  proceed:  First  transform  the  N UPPER  half 
spaces  to  N points  in  abc  space.  Construct  the  (bottom  part  of  the)  convex  hull  of  the 
points  in  OfNIogN)  time.  The  points  above  the  bottom  part  of  the  convex  hull 
correspond  to  “redundant’  half  spaces  and  can  be  discarded.  The  bottom  part  of  the 
hull  is  then  transformed  back  to  form  U,  the  intersection  of  the  UPPER  half  spaces. 
This  will  now  be  described  in  detail. 


, r j 

1 rU 


Figure  10:  A redundant  half  space. 


Assume  that  there  are  N UPPER  half  spaces.  Some  of  these  half  spaces 
contribute  to  the  intersection  U and  some  are  “redundant,"  such  as  half  space  M in 
Figure  10.  For  the  plane  there  are  two  simple  conditions  for  redundancy  of  a half 
plane.  In  three  dimensions  there  are  two  analogous  conditions  for  a half  space.  The 
first  condition  for  redundancy  of  a half  space  M is  that  plane  M be  above  the  point  P 
where  planes  I,  J,  and  K intersect,  as  in  Figure  10.  The  second  condition,  the 
"betweenness  of  slopes"  condition,  is  more  complicated  to  express.  The  purpose  of 
the  "betweenness  of  slopes"  condition  is  to  insure  that  a plane  above  the  point  P can’t 
drop  down  fas f enough  to  enter  the  region  below  planes  I,  J,  and  K.  But  this  is  insured 
by  requiring  that  plane  M remain  above  the  rays  r jj,  rj^,  and  r^j  which  define  the 
edges  of  the  three  faces  bounding  the  enclosed  region. 

Plane  M can  be  written  in  the  form  z - a*x  ♦ b*y  ♦ c.  (No  vertical  planes 
satisfying  the  two  conditions  are  possible,  so  this  form  will  always  work.)  Plane  I can 
be  written  z » aj*x  + b]*y  * cj.  Similarly  for  planes  J and  K.  The  three  rays  r j j,  r j^, 
and  rjq  can  be  expressed  as  vectors  originating  at  the  point  P.  The  cross  products  of 
normals  to  planes  I,  J,  and  K give  vectors  along  the  lines  determined  by  these  three 
vectors.  Thus, 
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Similarly  for  r and  r^j.  The  firs!  condition,  that  plane  M be  above  (or  on)  point  P - 
(P jj*P ylPj)  iS 

(1)  a*Py  + b*Pv  ♦ c i Pj. 


The  second  (betweenness  of  slopes)  condition  depends  only  on  the  slopes  a and  b,  not 
the  intercept  c.  Thus,  it  is  the  same  for  a plane  M which  passes  through  point  P as 


•i 

j 

one  which  is  strictly  above  P.  The  condition  that 
r^j  is  therefore 

M not  drop  below  vectors  rjj,  r or 

■i 

a * od,j  + 
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The  four  inequalities  in  (1)  and  (2)  have  an  equivalent  interpretation  in  the  transformed 
(abc) space. 


6.  Application  of  the  transform  to  the  three  dimensional  problem 


The  transform  used  in  three  dimensions  is  a straightforward  extension  of  the 
two  dimensional  transform.  Planes  transform  to  points  and  points  transform  to  planes. 
The  formulas  for  the  transform  are:* 


z - a*x  + b*y  + c -*  (a,b,c) 
(x,y,z)  -+  c - *x*a  ♦ -y*b  *■  z. 


*Dantzig[63]  uses  the  above  transform  in  the  context  of  linear  programming 
and  Huffman[77]  uses  an  almost  identical  transform  for  an  analysis  of  polyhedral 


scenes. 
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The  distance  between  a point  and  a plane  in  the  z coordinate  is  preserved  by  this 
transform.  And,  most  importantly,  the  sense  of  above/below-ness  (in  the  z coordinate) 
between  points  and  planes  is  preserved  also. 

Plane  M (z  - a*x  ♦ b*y  ♦ c)  of  Figure  10  transforms  to  the  point  (a,b,c)  in 
abc  space,  and  planes  I,  J,  and  K transform  to  the  points  Pj  ■ (aj.bj.cj),  Pj  - (aj,bj,r j), 
and  • (a^b^c^).  ^hc  inequality  (1)  (of  the  previous  section)  for  redundancy  of  half 
space  M (and  point  (a,b,c))  is  interpreted  in  abc  space  as  a half  space  whose  boundary 
plane  passes  through  all  three  points  Pj,  Pj,  and  P«.  Thus  condition  (1)  requires  that 
the  point  (a,b,c)  be  above  this  plane.  Inequalities  (2)  (of  the  previous  section)  define 
three  vertical  planes,  each  of  which  passes  through  two  of  the  three  points  Pj,  Pj,  and 
P^.  Thus,  the  set  of  halt  spaces  M’ which  half  spaces  I,  J,  and  K make  redundant  is 
represented  in  abc  space  by  the  set  of  all  points  (a,b,c)  which  are  directly  above  some 
point  in  (or  on)  the  triangle  determined  by  points  Pj,  Pj,  and  P|<.  When  conditions  (1) 
and  (2)  are  applied  over  all  N points  in  abc  space,  the  only  points  (half  spaces)  which 
are  not  made  redundant  by  some  points  (half  spaces)  I,  J,  and  K are  those  points  on 
the  bottom  part  of  the  conve*  hull  of  the  N points. 

How  long  does  it  lake  to  construct  the  (bottom  part  of  the)  conve*  hull  of  N 
points  in  three  dimensions’  Preparata  and  Hong(77]  describe  an  algorithm  to  construct 
the  (entire)  conve*  hull  in  OfNlogN)  time.  One  way  to  separate  the  bottom  of  the  hull 
from  the  top  requires  one  to  augment  Preparata  and  Hong’s  algorithm  by  maintaining 
with  each  face  of  the  convex  hull  a vector  perpendicular  to  the  face  which  points 
toward  the  inside  (as  opposed  to  the  outside)  of  the  hull.  The  bottom  faces  of  the  hull 
are  those  faces  whose  vectors  point  upward.  The  vectors  for  the  top  faces  point 
downward.  (Note  that  there  will  be  some  vertices  in  both  the  fop  and  bottom  par  Is  of 
the  hull.  These  are  the  vertices  which  bound  both  top  and  bottom  faces.)  The  bottom 
and  top  parts  of  the  convex  hull  are  therefore  separable  in  0(N)  time. 


Figure  11:  Transform  of  a convex  hull. 


To  find  the  intersection  of  the  UPPER  half  spaces,  the  bottom  part  of  the 
convex  hull  must  be  transformed  back  to  xyz  space.  The  0(N)  vertices  transform,  of 
course,  to  planes.  But  there’s  much  more  information  than  that  in  the  convex  hull.  For 
instance,  defining  each  face  of  the  hull  there  are  three  coplanar  vertices  I,  J,  and  K.l 
The  plane  which  these  vertices  define  is  transformed  to  a point  in  xyz  space.  This 
point  is  where  planes  I,  J,  and  K (of  xyz  space)  intersect.  (This  follows  from  Ihe  fact 
that  the  transform  preserves  incidence  between  points  and  planes.)  Also,  as  illustrated 


1 Assuming  that  the  vertices  are  in  general  position.  If  not,  then  there  may 
be  four  or  more  coplanar  vertices. 
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in  Figure  11,  if  faces  and  F^  of  the  (bottom  part  of  the)  convex  hull  share  an  edge 
VjVj,  then  in  xyz  space  faces  Vj  and  Vj  share  an  edge  F^F^.  In  fact,  even  the 
unbounded  faces  of  U in  xyz  space  can  be  obtained  from  the  transform.  These  fares 
of  U correspond  to  vertices  at  the  boundary  between  the  top  and  bottom  parts  of  the 
convex  hull  in  abc  space.  Thus,  little  computation  will  be  required  to  construct  U after 
the  transform  from  abc  space  since  all  the  faces,  vertices,  edges,  etc.  are  directly 
obtainable  from  the  transform.  Since  the  transform  costs  only  CXN)  time,  the  total  time 
to  construct  U is  dominated  by  the  time  to  construct  the  convex  hull  in  abc  space, 
which  is  (XNIogN)  time. 

It  has  just  been  shown  that  the  time  to  construct  U,  the  intersection  of  the 
UPPER  half  spaces  is  (XNIogN).  Thus,  U3,  L3,  U2,  (3.  Uj,  and  Lj  can  be  constructed  in 
O(NlogN)  time.  Once  the  Uj  and  Lj  are  constructed,  it  costs  only  0(N)  more  time  to 
intersect  U3  and  L3.  Thus,  the  total  time  to  intersect  N half  spaces  in  3 space  is 
CXNIogN). 

There  is  one  more  item  of  interest  in  the  three  dimensional  problem.  As 
with  the  two  dimensional  case,  an  intersection  of  half  spaces  algorithm  can  be  used  as 
a convex  hull  (of  points)  algorithm  also.  Simply  transform  all  N points  to  UPPER 
(LOWER)  half  spaces,  intersect  the  half  spaces,  and  transform  back  to  obtain  the 
bottom  (top)  part  of  the  convex  hull.  Merging  the  top  and  bottom  parts  of  the  hull  is, 
again,  trivial  since  they  will  share  several  vertices  at  the  boundary  of  the  top  and 
bottom  parts. 


7.  Intersecting  halt  spaces  in  four  or  more  dimensions 

Suppose  that  in  K dimensions  the  convex  hull  of  N points  can  be  constructed 
in  H(N,K)  time.  Then  the  intersection  of  N half  spaces  can  be  constructed  in  0(K*H(N,K)) 
time.  The  problem  is  solved  by  a straightforward  extension  of  the  three  dimensional 
solution.  The  Uj  and  Lj  are  first  constructed  in  (XK*H(N,K))  time  and  then  U«  and  L« 
are  intersected  with  the  aid  of  the  other  Uj  and  Lj  as  in  the  three  dimensional  case.  To 
construct  the  Uj  and  Lj,  one  transforms  the  UPPERj  (LOWERj)  half  spaces  to  points, 
constructs  the  convex  hull  of  the  points  in  H(N,K)  time,  partitions  the  top  from  the 
bottom  parts  of  the  hull,  and  transforms  the  bottom  (top)  part  of  the  hull  back. 

Before  proving  that  this  actually  works  in  K dimensions  some  terminology 
will  first  have  to  be  introduced.  Lef  the  UPPER  set  of  half  spaces,  refer,  by  default,  to 
the  UPPERj  set  of  half  spaces.  That  is,  "up"  will  be  measured  in  the  xj<  coordinate. 
Also,  let  a *j-face"  of  a K dimensional  polytope  by  denoted  as  follows:  a vertex  is  a 0- 
face,  an  edge  (a  line  segment)  is  a 1-face,  etc. 


Theorem:  Let  S be  a set  of  N UPPER,  K - dimensional  half  spaces  Hj  defined  by 
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and  let  T be  the  convey  polytope  formed  by  the  intersection  of  these  N half  spaces.  Jf 
the  transform  of  half  space  Wj  r S is 

Hj  -*  Aj  - <i»j i .ai2—aik)» 

then  the  nonredundant  half  spaces  of  S correspond  to  the  points  on  the  bottom  half  of 
the  convex  hull  of  the  points  Aj.  Furthermore,  the  j-faces  of  the  convex  hull 
correspond  to  the  K-j-l  faces  of  T. 

Proof:  The  proof  is  divided  into  five  parts: 

(1)  For  each  redundant  half  space  Hj  there  is  a set  Gj  c (S-Hj)  of  half  spaces  such 
that  [ INTERSECT(x'Gj)  x ] c Hj  and  |Gt|  - K (if  the  half  spaces  are  in  general 
position). 

(2)  The  algebraic  description  of  the  redundancy  of  H|  with  respect  to  Gj. 

(3)  The  interpretation  of  the  algebraic  description  in  the  transformed  space. 

(4)  The  interpretation  in  (3)  implies  that  the  nonredundant  half  spaces  correspond  to 
the  points  on  the  bottom  part  of  the  convex  hull  of  the  points  Aj. 

(5)  The  j-faces  of  the  convex  hull  correspond  to  the  K-j-l  faces  of  T. 


PART  (1):  T is  a convex  polytope  because  it  is  an  intersection  of  half  spaces.  If  a 
halfspace  Hj  is  redundant  then  its  boundary  plane  lies  completely  above  T.  Let  Vj  be 
the  point  of  T which  is  closest  to  the  boundary  plane  of  Hj.  (If  the  closest  point  to  the 
boundary  of  Hj  is  not  unique,  then  let  Vj  be  any  vertex  in  the  set  of  closest  points.) 
Since  T is  a convex  polytope,  Vj  is  (or  can  be  chosen  to  be)  a vertex  of  T.  Let  Gj  c s 
be  the  set  of  half  spaces  of  S whose  boundaries  meet  at  vertex  Vj.  If  the  half  spaces 
of  S are  in  general  position  thch  the  set  G|  will  contain  exactly  K half  spaces.  Since  T 
is  convex,  one  can  travel  from  Vj  along  T in  any  direction  and  the  distance  to  the 
boundary  of  Hj  will  be  nondecreasing.  It  follows  that  [ INTERSECT(xcGj)  x ] c Hj.  That 
is,  Hj  is  redundant  with  respect  to  the  half  spaces  in  Gj.  As  in  two  and  three 
dimensions,  redundancy  can  be  expressed  by  two  sets  of  conditions:  (1)  the  boundary 
of  Hj  is  above  vertex  Vj,  and  (2)  the  “betweenness  of  slopes"  condition.  These 
conditions  will  now  be  shown  algebraically. 

PART  (2):  Let  Hj  be  Ihe  half  space 

«»k  4 2-  j*Xj  S xk, 

TK  j-l.K-1  ,J  1 * 

and  let  the  half  spaces  in  Gj  be,  without  loss  in  generality,  the  half  spaces  1 through  K: 

»ik  4 aL  aii**i  * xk>  i”1— K- 
j-UM  1 J 

The  point  Vj  - (Vjj,Vj2,  • • • Vfk)  closest  to  the  boundary  plane  of  Hj  is  defined  as  the 
solution  to  the  set  of  equations 
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Letting 


A 


aik  * ^ " ^tk' 

J"  1 >K-I 


this  can  be  rewritten 


A » V,’  - B, 

where  Vj  is  the  transpose  ol  Vj. 

Now  the  two  conditions  lor  redundancy  of  Hj  can  be  described.  The 
condition  that  the  boundary  ol  H|  is  abovo  point  Vj  is 

(1>  «tk  ♦ . JE  VVH  2 Vtk* 

j-l.K-l 

The  "betweenness  ol  slopes"  conditions  are 

(2)  X MjjiCi  2 0,  i-l,..K, 

j-l.K  1 ' 

where  M(j  - (-1)K  * cofjj(A)  * SGNfcofjK(A))  and 
Cj  - ajj  lor  j<K-l  and  - -1. 

Equation  (I)  is  straightforward  and  believable,  but  (2)  requires  some 
explanation.  Describing  (?),  (he  "betweenness  of  slopes*  condition,  is  difficult  to  do  by 
directly  comparing  slopes  of  the  K flats  of  Gj.  Instead,  the  alternative  developed  lor 
the  three  dimensional  case  will  now  be  applied  to  the  K dimensional  problem.  In  three 
dimensions  one  constructs  the  three  rays  rjj,  r^,  and  r^  determined  by  the  vertex  P 
where  the  three  planes  meet  and  the  intersection  ol  pairs  of  planes  I and  J,  I and  K, 
and  J and  K,  respectively.  Condition  (2)  is  satisfied  if  the  redundant  plane  stays  above 
these  rays.  In  K dimensions  one  constructs  the  K rays  determined  by  the  vertex  Vj 
and  the  K sets  of  K-l  elements  of  G|.  The  notation  used  for  the  K dimensional  case 
will  be  slightly  different  than  lor  the  three  dimensional  case.  Rather  than  denoting  a 
ray  by  the  set  of  K-l  flats  which  determine  it,  the  ray  will  instead  be  denoted  by  the 
one  flat  which  Is  not  used  to  determine  it.  If  (<*j]>otj2-  . . . is  a vector  parallel  to 
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ray  i (the  ray  which  docs  not  include  flat  i),  then  ray  i may  be  written  in  parametric 
form  as 


(*1»*2*  • • • xk^  " ^tl>^t2‘  • • • ^tk^  u # (uHle^i2,  • • • rtik^ 

Now  It  must  be  shown  how  to  construct  such  a vector  (oCj  j >c^-i2*  • • • °^ik^  ' 's  also 

the  intersection  of  the  other  K-l  half  spaces  of  Gj.  Thus, 


By  combining  these  two  equations  one  obtains 


The  problem  is  now  to  find  a general  solution  to  this  system  of  K-l  equations  in  K 
unknowns.  Recall  the  following  property  of  cofactors: 

X ajj*cofm:(A)  - det(A)  if  i-m 
i-l.K  - 0 if  i K m. 

It  follows  that  the  general  solution  is 


. . . ociK)  - constant  * (cofjj(A),cofj2(A), . . . cofjk(A)). 


I 
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The  constant  is  chosen  as  (-1)*  • SGN(cofj^(A))  so  that  the  vectors  always  point 
downward.  If  can  now  be  seen  that  the  inequality  expressed  in  condition  (2)  is  simply 
a generalization  ot  the  formulas  lor  the  three  dimensional  case. 

PART  (3):  In  the  transformed  space,  half  space  Hj  becomes  a point  Hj  - (atl>',t2’  • 
. . aij^)  and  the  elements  of  Gj  transform  to  points  similarly.  Vertex  Vj  transforms  to  a 
flat  V{  defined  by 

V|k  ‘ j.§.i  V|i‘'i  * 'k- 

Since  incidence  is  preserved  by  the  transform,  flat  Vj  passes  through  all  of  the  points 
Gj  and  is,  in  fact,  determined  by  them.  Also,  since  above/below-ness  in  the 
coordinate  is  preserved,  point  Hj  is  above  flat  Vj  in  the  transform  space.  Using  the 
formula  for  flat  Vj,  this  can  be  expressed  as 

v,k  ■ i-S-t Va,i  s ,|k' 

which  is  equivalent  to  condition  (1)  for  redundancy  of  Hj!  Condition  (2)  defines  K 
"vertical"  half  spaces  which  further  constrain  point  Hj  to  be  not  only  above  flat  Vj  but 
also  directly  above  some  point  in  the  (K- 1 dimensional)  simplex  of  the  points  Gj.  (The 
K half  spaces  are  described  as  "vertical"  because  their  formulas  do  not  include  the 
term  a^  and  thus  the  K coordinate  is  free.) 

PART  (4):  In  Parts  (!)  through  (3)  attention  has  focused  only  on  the  conditions  for 
redundancy  of  a particular  half  space  Hj.  Given  the  interpretation  of  part  (3)  above, 
what  can  now  be  said  over  all  N of  the  half  spaces  of  S?  It  has  been  shown  that  if 
half  space  Hj  is  redundant  with  respect  to  a set  of  K half  spaces  Gj,  then  in  the 
transform  space  point  Hj  is  directly  above  some  point  of  the  (K-l  dimensional)  simplex 
of  the  points  Gj,  and  conversely.  Thus,  a hall  space  Hj  ( S is  redundant  iff  there  exists 
a set  Gj  c S-Hj  of  K halt  spaces  such  that  in  the  transform  space  point  Hj  is  directly 
above  a point  of  the  (K-l  dimensional)  simplex  of  the  points  of  Gj!  But  this  eliminates 
all  points  except  those  on  the  bottom  of  the  convex  hull  of  the  N points  Hj.  This 
proves  the  main  part  of  the  theorem.  Now  it  must  be  shown  that  the  convex  hull 
produces  not  only  the  nonredundant  half  spaces  but  also  the  edges,  faces,  etc.  of  T 
also. 

PART  (5):  To  show  that  a j-face  of  the  bottom  part  of  the  convex  hull  of  the 
points  Hj  corresponds  to  a K-j-l  face  of  T,  a generalized  version  of  the  transform  must 
first  be  described.  The  general  transform  maps  a j-space  to  a K-j-1  space. 

General  Transform:  The  transform  of  a j dimensional  subspace  of  K 
space  is  the  set  of  all  flats  which  contain  it.  This  is  a K-j-1  dimensional 
subspace  of  flats.  However,  since  each  flat  can  be  readily  represented  as 
a point,  the  transform  of  the  j dimensional  subspace  is  represented  as  a 
K-j-1  dimensional  subspace  of  points. 


Theorem:  The  general  transform  preserves  incidence.  In  other  words,  if  a jj 
dimensional  subspace  ot  K space  is  a subspace  of  a 12  dimensional  subspace,  then  the 
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transform  of  the  jj  dimensional  subspace  is  a subspace  of  the  transform  of  the  jj 
dimensional  subspace 

Proof:  Since  the  12  dimensional  subspace  can  be  interpreted  as  an  intersection  of  K-j2 
flats,  it  will  be  sufficient  to  show  that  if  a K-j  dimensional  subspace  is  a subspace  of  a 
flat,  then  the  transform  of  the  flat  is  a (zero  dimensional)  subspace  of  the  transform  of 
the  K-j-space.  Let  the  K-j-space  he  represented  as  an  intersection  of  j flats 


It  is  easy  to  see  that  this  j-l  space  contains  the  points 

<aj i.»j2*  • • • •i,k-l>"ik)* 
which  are  the  transforms  of  the  j flats. 

Since  the  general  transform  preserves  incidence  and  j-spaces  transform  to  K-j-1 
spaces,  it  follows  that  j-faccs  of  the  convex  hull  transform  to  K-j-1  faces  of  T.  QED. 


As  with  the  two  and  three  dimensional  cases,  an  algorithm  for  intersecting  N 
half  spaces  in  K dimensions  can  be  used  to  construct  the  convex  hull  of  N points  in  K 
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dimension*.  The  procedure  i*  a straightforward  extension  of  the  three  dimensional 
case:  Transform  all  N points  to  N UPPER  (LOWER)  half  spaces  in  K space.  Intersect  the 
half  spaces,  then  transform  back  to  the  LOWER  (UPPER)  part  of  the  convex  hull. 
Joining  the  UPPER  and  LOWER  parts  of  the  hull  should  be  easy  because  of  the  vertices 
which  they  have  in  common  on  their  boundaries. 


8.  Fast  expected  time  algorithms 

There  have  been  recent  advances  in  algorithms  for  constructing  the  convex 
hull  of  a set  of  points  in  fast  expected  time.  Since  the  time  to  intersect  N half  spaces 
is  dominated  by  the  time  to  construct  the  convex  hull  of  N points,  these  fast  convex 
hull  algorithms  lead  to  fast  intersection  of  half  space  algorithms.  Bentley  and 
Shamos[78]  have  shown  that  when  the  expected  number  of  points  on  the  convex  hull 
is  O(NP)  for  some  p < 1,  the  convex  hull  of  N points  in  two  and  three  dimensions  can 
be  constructed  in  0(N)  expected  time,  while  maintaining  an  (XNIogN)  worst  case  lime. 
Thus,  N half  planes  or  N half  3-spaces  (under  the  above  conditions)  can  be  intersected 
in  0(N)  expected  time  and  QfNIogN)  worst  case  time.  ^ In  four  or  more  dimensions  the 
worst  case  time  is  at  least  O(N^).  However,  it  is  anticipated  that  one  may  still 
construct  fast  expected  lime  algorithms. 


9.  Applications 

(A)  Linear  programming. 

The  feasible  region  of  a linear  programming  problem  in  K variables  and  N 
constraints  is  an  intersection  of  N K-dimensional  half  spaces.  For  large  K,  the  time 
0(K*H(N,K))  to  construct  this  may  not  be  an  improvement  over  the  simplex  method.  But 
for  K S 3,  the  feasible  region  can  be  constructed  in  CXNIogN)  time  worst  case  and 
(under  the  conditions  described  in  Section  8)  CXN)  expected  time.  Thus,  the  total  time 
to  solve  the  linear  programming  problem  is  only  0(NlogN)  worst  case  and  CXN) 
expected  time.  Shamos  and  Hoey[76]  have  previously  shown  that  a linear 
programming  problem  in  two  variables  (two  dimensions)  can  be  solved  in  0(NlogN)  time. 
Bentley  and  Shamos[73]  later  showed  the  0(N)  expected  time  for  the  two  dimensional 
case. 

For  K > 3 dimensions  a different  approach  is  recommended.  To  solve  a 
linear  programming  problem  it  is  not  really  necessary  to  construct  the  entire  feasible 


1 Bentley  and  Shamos[78]  showed  this  result  for  the  two  dimensional  case 
using  a duality  result  of  Ziezold[70J. 
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region.  Certainly  the  simplex  method  avoids  doing  any  such  thing.  Dantzig[63]  applies 
the  duality  transform  to  a standard  linear  programming  problem  of  K equations  in  N 
. variables.  (The  interpretation  of  a linear  programming  problem  as  an  intersection  of  N 

half  spaces  in  K dimensions  is  shown  equivalent  to  this  standard  form  by  the  duality 
theorem  of  linear  programming.)  The  objective  function  is  transformed  to  a vertical  line 
in  K space  and  the  N variables  are  each  transformed  to  points  in  K space.  The 
optimization  problem  is  transformed  to  the  problem  ol  determining  where  the  vertical 
line  intersects  the  convex  hull  of  the  N points.  The  simple*  method  provides  one  way 
to  determine  this  intersection  without  constructing  the  convex  hull  of  the  N points.  It 
remains  to  be  determined  whether  or  not  modern  algorithm  techniques  can  produce  an 
improvement  over  the  simplex  algorithm  in  the  K dimensional  case. 

(B)  Intersection  of  convex  polyhedra 

A convex  polyhedron  of  N faces  can  be  described  as  an  intersection  of  the  N 
half  spaces  which  determine  the  faces  (and  which  include  the  polyhedron).  Thus,  an 
intersection  of  two  convex  polyhedra  of  N faces  can  be  solved  as  an  intersection  ol 
2*N  half  spaces,  which  can  be  done  in  OfNIogN)  time.  (See  Muller  and  Preparata[77] 
for  a different  OOMIogM)  time  algorithm  for  intersection  of  convex  polyhedra.)  But  it 
remains  an  open  question  whether  or  not  the  algorithm  can  be  improved  to  run  in  0(N) 
time.  (It  certainly  can  in  two  dimensions.) 

' 
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10.  Conclusion 

The  use  of  a geometric  transform  has  been  shown  to  be  very  useful  for 
construction  of  a fast  algorithm  for  intersection  of  half  spaces.  This  algorithm,  in  turn, 
has  been  used  in  algorithms  for  a fast  expected  time  algorithm  for  intersection  of  half 
spaces,  intersection  ol  convex  polyhedra,  and  linear  programming  in  three  variables. 
The  geometric  transform  is  applicable  to  many  more  problems.  The  union  of  a set  of 
(the  interiors  of)  N arbitrary  planar  circles  can  be  constructed  by  transforming  to  a set 
of  N half  3 - spaces  which  are  then  intersected  in  O(NlogN)  time.  The  Euclidean 
diameter  of  a set  of  N points  in  3 - space  can  be  determined  in  OfNIogN)  time  through 
use  of  a transform  very  similar  to  the  one  in  this  paper.  The  nearest  and  farthest 
point  Voronoi  diagrams  ol  N points  on  a sphere  or  on  a Euclidean  plane  can  be 
constructed  in  0(NlogN>  time  through  use  ol  inversion  and  the  methods  in  this  paper. 
Also,  several  two  dimensional  problems  with  lines  have  been  solved  quickly  through 
use  of  this  transform  <Brown[77]).  In  addition  to  the  duality  transform  described  in 
this  paper,  other  transforms  such  as  inversion,  rotation,  and  various  kinds  ot  projection 
(gnomonic,  orthographic,  ster cographic)  have  been  successfully  used  to  construct  fast 
geometric  algorithms.  It  is  anticipated  that  the  use  of  geometric  transforms  will 
become  a standard  tool  in  the  construction  of  geometric  algorithms. 
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